Macrodimers: ultralong range Rydberg molecules 
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We study long range interactions between two Rydberg atoms and predict the existence of ultra- 
long range Rydberg dimers with equilibrium distances of many thousand Bohr radii. We calculate 
the dispersion coefficients Cs, Ce and Cs for two rubidium atoms in the same excited level np, and 
find that they scale like n*, n^^ and n^^ , respectively. We show that for certain molecular symme- 
tries, these coefficients lead to long range potential wells that can support molecular bound levels. 
Such macrodimers would be very sensitive to their environment, and could probe weak interactions. 
We suggest experiments to detect these macrodimers. 

PACS numbers: 34.20.Cf, 32.80.Rm, 32.80.Pj 
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New techniques used for cooling and trapping of atoms 
and molecules and which led to the realization of 
atomic Bose-Einstein condensation have also been 
applied to experiments with ultracold plasmas and 
ultracold Rydberg atoms [||. The exaggerated proper- 
ties of Rydberg atoms provide a fertile ground for new 
physics. For example, transport properties of ultracold 
gases doped with ions were recently explored and ex- 
tended to cold Rydberg samples ||^, while entangled 
states relevant for quantum computing can also be pro- 
duced with ultracold Rydberg atoms 0. Finally, the 
creation of "trilobite" Rydberg molecules was proposed 
1^ , where one atom of the dimer remains in its ground 
state while the second one is excited to a Rydberg state. 

In this paper, we explore the interactions between two 
Rydberg atoms. We show that long range wells support- 
ing several bound levels exist for certain molecular sym- 
metries. We explore the sensitivity of these wells to the 
particular asymptotic form of the potential long range ex- 
pansion and show that their existence is robust. We also 
estimate the effect of retardation as well as the validity 
of the Born-Oppenheimer approximation. We give nu- 
merical examples for the case of rubidium (Rb) , and dis- 
cuss experimental schemes to detect these macrodimers. 
These molecules could lead to measurements of very weak 
interactions, such as vacuum fluctuations, and provide a 
unique tool to study quenching in ultracold collisions. 

We consider two atoms each excited by one photon 
from their ground state into the same Rydberg state np 
1^, where n is the principal quantum number. For Rb, 
this corresponds to the 5s — > np transition jl^. At large 
separation R, the potential energy between two atoms 
can be expanded in powers of 1 /R |0 . For two identical 
atoms in the same np state, it takes the form [|2| 
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where the dispersion coefficients C5, Cg, and Cg depend 
on n. For the np — np asymptote of the homonuclear 
dimers, we have in total six pairs of degenerate molecular 
states with identical coefficients for each pair: 
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FIG. 1: Comparison of the potential curves for the three 
pairs of states that sustain a well for Rb(rj = 20). Inside the 
shaded area we have R < -Rlr, and Eq. (Ilh is not adequate. 



^A„, ^Ilg and ^n„, ^Ilg and ^n„, and ^S", and two 
pairs of and Note that the degeneracy of 

any pair of molecular states is lifted by the exponentially 
decaying exchange interactions, which can be neglected 
at large enough R. 

The numerical values of C5, Cg and Cs |Q were cal- 
culated using the expressions of Marinescu [Q. The 
sums over the electronic states were evaluated directly, 
and using the rt-scaling of the dipole (oc r?\ quadrupole 
(cx n^), and octopole (oc n^) matrix elements and the 
energy differences (oc n~^) involved, we obtained the fol- 
lowing scaling laws: C5 oc n®, Cg oc n^^, and Cg oc n}'° . 
The magnitude and sign of the dispersion coefficients de- 
pend on the molecular symmetry considered, and it is 
possible to obtain a long range potential well with an 
attractive long range R~^ contribution and a repulsive 
shorter range R~^ or R~^ contribution (see Fig. |^). For 
the np — np asymptote, we found that three pairs of de- 
generate molecular states give long range potential wells: 
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-^S^j , and one 
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For the system to be adequately described by Eq. (|^), 
the exchange energy must be negligible. To estimate 
the region of vahdity of Eq. (§), we use the Le Roy 
radius i?LR as measure of the electron wavefunc- 
tion overlap between the two atoms: it is given by 
i?LR = 2 (y^fl^ -)- (r2)^/2^, where (r^)^/^ is the rms 
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FIG. 2: In (a), comparison of De as a function of Re for the 
three pairs supporting long range wells, for various n. In (b), 
iie is compared to _Rlr for the same states. The shaded area 
is defined by Re < RhR- 

position of the electron of atom A [B) If i? < i?LR, 
exchange and charge-overlap interactions become impor- 
tant and Eq. (|]) is not adequate to describe the system. 

In Fig. |lj we illustrate the wells for the three pairs 
of degenerate states ^Hg-^n^, ^T^'-^Y.-, and 
of Rb with n = 20. The Le Roy radius i?LR is equal to 
1902 ao, where oq is a Bohr radius. The well of the pair of 
states ^E+-^I]+ is very deep (potential depth — 4.5 x 
10^ cm~^), but is located at a much shorter distance than 
the Le Roy radius (equilibrium distance Re = 240 oq): 
the interactions for these states are not well described by 
Eq. (|l|), as opposed to the two other pairs. The well of the 
pair ^Ilg-^IIti has a depth = 3.74 x 10~^ cm"'^ and an 
equilibrium distance Re — 2509 cq. By comparison, the 
well of the pair '^E~-^S~ is much shallower and farther 
away, with De = 6.88 x 10~^ cm~^ and Re — 4956 aq. 
In Fig. H, we compare ZJg, Re and i?LR for the three 
pairs for various values of n, and find the same general 
behavior. For the remainder of this letter, we focus our 
attention on the deeper wells described by Eq. (|l|), i.e., 
the "'^ng-'^n„ pair, since the much shallower wells of the 
pair may prove more difficult to detect. 

In Fig. m we illustrate the scaling of the dispersion 
coefficients, equilibrium distance, and well depth of the 
^rtg-'^IIi, pair as a function of n. In atomic units, the 
coefficients are given approximately by C5 ^ 3 n^, 
Ce 0.7 n", and Cg 50 n^^ (see Fig. |(a)). Ne- 
glecting Cs in Eq. (p, and setting the derivative of V{R) 
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FIG. 3: Scaling of (a) C5, Ce, Cg for the ^n<,-^n„ pair, and 
(b) Re (left scale) and De (right scale), as a function of n. 



to zero, we find that the equilibrium distance scales as 
Re 0.3 ao, in good agreement with the numerical 
values shown in Fig. ^(b). The scaling of De is 10^n~^ 
cm^^, as expected (see Fig. ||(c)). Note that other pairs 
may have different Re and De scaling, depending on the 
relative magnitude of the dispersion terms. The wells for 
the ^Hg-^Uu pair, although shallow, support many vibra- 
tional iDOund levels. In Table |, we list the two lowest and 
highest levels found for n = 20, 40 and 70 which sup- 
port 143, 125 and 107 bound levels, respectively. While 
the wells for high n are much shallower than those with 
smaller n (e.g., see n = 20 and 70 in Table |), their larger 
extension leads to denser energy levels, and hence they 
also support a large number of levels. 

When using the expression for V{R), the Born- 
Oppenheimer approximation is assumed valid. To ver- 
ify that it is the case, we compare the vibrational period 
Tvib (v) of a given bound level v with the typical time for 
the electron motion Tg = h/\En\ = {4:Tr / acjn^ ao where 
En is the electron binding energy. If Te/rvib <S 1, the Ry- 
dberg electrons of both atoms easily follow the motion of 
the ions, and the Born-Oppenheimer curve V{R) needs 
no further diabatic corrections. Taking the levels in Ta- 
ble|| as examples, for n = 20, we find Te/rvib ~ 1-5 x 10~^ 
for the lowest level v = 0, and 2.2 x 10^^^ for the highest 
level V = 142: the electrons follow adiabatically the ions 
for all vibrational levels. For n = 70, the ratio becomes 
6.8 X 10'^ for = and 4.5 x 10"" for v = 106. Natu- 
rally, the Born-Oppenheimer approximation gets better 
for the higher vibrational levels, since they are more ex- 
tended and therefore take more time to make a full os- 
cillation. However, the spontaneous decay of one of the 
excited Rydberg atoms will limit the lifetime of these long 
range molecules and prevent the existence of the upper 
lying vibrational bound levels. Although the lifetime of 
the excited atoms is long, scaling as Tat Ton^, where 
To ^ 1.4 ns it is much shorter than the vibrational 
period for high v. E.g., for n = 20, we find Tat ~ 11-2 
/iS and Tvib = 0.083 ^s and 5.58 ms for w = and 142, 
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TABLE I; Sample of vibrational bound levels for the ^Ilg- 
pair. The top corresponds to n = 20, with Re = 2509 ao, 
= 3.740 X 10"^ cm"\ the middle to n = 40, with = 
20141 ao, De = 4.060x10"" cm"\ and the bottom to n = 70, 
with Re = 103 380 ao, De = 1.042 x 10"^ cm"^ 
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-3.70 [- 


^2] 


4.01 [-4] 


2 454 


2 571 


8.3 


[-8] 
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-3.62 [- 


-2] 


1.19 [-3] 


2417 


2 623 


8.4 


[-8] 


141 


-3.46 [- 


^8] 


3.74 [-2] 


2142 


49 859 


1.2 


[-3] 


142 


-3.89 [- 
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3.74 [-2] 


2142 


77 220 


5.6 


[-3] 
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19 651 


20 701 
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[-6] 
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-3.91 


[-4] 


1.50 
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19 331 


21163 
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[-6] 
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-3.27 [- 
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17111 


431 262 
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-1.91 [- 
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98 762 
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[-4] 


105 


-1.02 [- 


-11] 


1.04 


-5] 


87 306 


2 211658 




3.1 
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-3.57 [- 


-13] 


1.04 


"5] 


87 306 


4 337482 




32.8 



respectively, and for n = 70, Tat ^ 480 /iS and Tvib = 220 
fis and 32.8 s for u = and 106, respectively. In ef- 
fect, the upper vibrational levels can be considered as 
quasi-continuum states since only a fraction of an entire 
oscillation will take place before de-excitation. 

Another effect limiting the existence of Rydberg long 
range molecules is auto-ionization. When two excited 
atoms interact, one atom can decay to a lower excited 
state while the second atom is ionized, the free elec- 
tron picking up most of the kinetic energy. However, 
if the separation between the Rydberg atoms is larger 
than RhR, there is little overlap between their electronic 
clouds, and one expects the auto-ionization probability 
to be small (the atoms interact via their electric dipoles 
at such distances). Note also that for alkali dimers in 
general [p^ , and Rb in particular, avoided crossings with 
the ionic curves (here correlated to Rb++Rb^) perturb 
the lowest asymptotes such as bp — bp. However, for the 
case where the Le Roy radius condition is satisfied, i.e., 
for n > 20 (see Fig. ^(b)), the np — np asymptotes are 
much higher than the ionic curves, and no perturbation 
from avoided crossings will occur. 

The density of electronic states in the Rydberg region 
is very large, and avoided crossings from potential curves 
correlated to other asymptotes with the same symmetry 
could perturb the shallow wells. However, the closest 
asymptotes are generally separated by several times the 
potential depth Dg, and no avoided crossing in the well 
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FIG. 4: Nearest asymptotes to 40p — 40p for the ^Ug-^JJu 
symmetries. Inset: convergence of Eq.(^ for n = 40. The 
shaded area gives an estimate of the truncation error. 



region occur . Fig. ^ illustrates the well for 40p — 40p 
and the two nearest curves correlated to 36s — 47p and 
42s — 36/, for the ^Hg-^Hu symmetries. The separa- 
tion between them is much larger than and no cross- 
ings occur. Although the various asymptotes become 
denser as n increases, the potentials become very shal- 
low (Z?e oc and generally no crossing is found. Note 
that avoided crossings from a lower asymptote repulsive 
curve and a higher asymptote attractive curve (with the 
same symmetry) could, in fact, produce deep wells. 

For these extremely extended states, retardation ef- 
fects may become important. Again, to estimate the 
importance of the photon time-of-flight, we compare the 
electronic time Tg with the time it takes a photon to cover 
the distance between the two centers: Tph ~ R/c. For the 
highest levels, the relevant distance is the outer turning 
point i?2, but for the deeper levels, a good estimate is ob- 
tained by using the equilibrium distance Rg. Using the 
scaling R^ ^ 0.3 n^, we have Tph/rg ^ (0.3a/47r)n ^ 1 
for n < 100. For the highest vibrational levels v in Ta- 
ble I we get Tph/Te - 0.06,0.07, and 0.11, for i?2 of 
V — 140,141, and 142, respectively (for n — 20), and 
Tph/Te ~ 0.19, 0.26, and 0.51 for i?2 of ^; = 104, 105, and 
106, respectively (for n = 70). Although retardation ef- 
fects are not important for the lower vibrational levels 
(with i?2 slightly larger than Re), even for large n, one 
needs to consider them for the higher vibrational levels. 
The sensitivity of the interaction potential to retardation 
effects could actually be used to detect these effects with 
high accuracy using high precision spectroscopy. 

We included three terms in the asymptotic 1 /R expan- 
sion (P , but one has to verify the effect of the truncation 
on the well properties. The standard procedure to esti- 
mate the maximum error due to the truncation of higher 
terms is to include half of the last term here, C's- In 
Fig. 1^, we illustrate the effect of such variation for n = 40. 
The well becomes deeper, but is still located outside the 
Le Roy radius, hence the existence of long range wells is 
robust to the truncation of the 1/R expansion. 
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To detect these Rydberg macrodimers, one could pho- 
toassociate them from the ground state (5s for Rb) to 
the target molecular symmetry, using appropriate pho- 
ton wavelength and polarization. After the photoassoci- 
ation pulse pO[ |, one can field-ionize the system with an 
electric pulse and detect the ions. The existence of the 
molecules should lead to small peaks red-detuned from 
the atomic asymptote. Because the spacing of the bound 
levels is small, individual levels will be hard to detect (a 
laser width smaller than the spacing A is needed: for 
n = 20, A - 10^3 cm-i ~ 300 MHz for the deepest 
levels), and the most probable signature would be an 
overall red-detuned peak. The wells are shallow, and 
only those deep enough could be detected. To enhance 
the sensitivity of the detection scheme, one could also 
use radio-frequency to excite a bound level to a contin- 
uum state corresponding to another asymptote (e.g., to 
{n + l)s — np). The detection of any {n + l)s atoms (e.g., 
by ramp-field ionization) would be attributed to the ex- 
istence of a long range bound molecular state, since the 
RF field would not be resonant with free np atoms. 

We have computed long range interactions between two 
Rb atoms in the same np state, and found that shallow 
long range wells supporting several vibrational levels ex- 
ist for certain molecular symmetries. Although specific 
calculations were performed for Rb and np — np asymp- 
totes, the existence of these macrodimers is general: one 
can expect them for various asymptotes niii — 71.2^2 and 
for all alkali atoms. Furthermore, avoided crossing be- 
tween curves of the same symmetry with different asymp- 
totes could also provide deep long range molecular wells, 
in a manner similar to the long range wells observed in 
many alkali dimers . The macrodimers are extremely 
sensitive to their environment, and as such, they could be 
used as probes for extremely weak interactions, e.g., to 
measure retardation effects and vacuum fluctuations, or 
any weak electromagnetic interaction. Also, due to their 
very small rovibrational energy splittings, macrodimers 
would provide a unique tool to study quenching in ul- 
tracold collisions, since they would remain trapped (as 
opposed to usual dimers [^). Finally, the detection of 
such exotic molecules in itself would be a considerable 
achievement. 

The authors thank A. Dalgarno, V. Kharchenko, P.L. 
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ITR-0082913 from the National Science Foundation. 
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